Odd-Parity Triplet Pair Induced by Hund's Rule Coupling 
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We discuss microscopic aspects of odd-parity triplet pair in orbital degenerate systems. From the 
concept of off-diagonal long-range order, a pair state is unambiguously defined as the eigenstate 
with the maximum eigenvalue of pair correlation function. Performing this scheme by a numerical 
technique, we clarify that the odd-parity triplet pair occurs as an out-of-phase combination of local 
triplets induced by Hund's rule coupling for the lattice including two sites in the unit cell. 

PACS numbers: 74.20.Rp, 74.20.Mn, 74.70.Tx, 74.70.Pq 
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Recently triplet superconductivity has been found in 
several compounds with strong electron correlation. A 
pioneering material should be UPts jlj, which has been 
considered to have spin-triplet pairing from experimen- 
tal evidence for multi superconducting phases similar to 
the phase diagram of "^He. The odd-parity triplet pair 
has been eventually confirmed by the experimental fact 
that the Knight shift does not change through the su- 
perconducting transition temperature 2]. In other /- 
electron compounds, triplet superconductivity has been 
also suggested. For instance, coexistence of superconduc- 
tivity and ferromagnetism has been discovered in UGe2 
[3 and URhGe Q . It is naively believed that spin triplet 
pair appears in the ferromagnetic (FM) phase. Recently 
UNi2Al3 has been also considered as a triplet supercon- 
ductor from NMR measurements 0]. 

Among transition metal oxides, Sr2Ru04 has attracted 
much attention since its discovery in 1994 It is 

confirmed that this material has spin triplet pair from 
NMR measurements 0], although it is isostructural to 
La2Cu04 which is the mother compound of high-Tc 
cuprates with singlet c?-wave superconductivity. In the d- 
electron system such as ZrZn2, it has been also reported 
that superconductivity appears in the FM phase |3| and 
spin-triplet pairing is expected to occur in this material. 

Regarding the mechanism of triplet superconductivity, 
FM spin fluctuation was considered to mediate the triplet 
pair also in strongly correlated materials by analogy with 
superfluid '^He. However, paramagnons are not always 
dominant in the spin fluctuation spectrum of those ma- 
terials. In fact, in Sr2Ru04, significant enhancement of 
the incommensurate antiferromagnetic spin fluctuation 
has been observed in neutron scattering experiments 
Namely, in contrast to the naive expectation, param- 
agnons do not always play a central role in the occur- 
rence of triplet superconductivity in the solid state. It is 
still a puzzling and challenging problem to clarify a key 
issue to determine the Cooper-pair symmetry. 

One possible scenario is based on the Hund's rule cou- 
pling, which stabilizes the local triplet pair composed of a 
couple of electrons between different orbitals. Then, the 
triplet superconductivity is naively expected to occur. 



but in such a local picture, it is questionable whether the 
anisotropic Cooper-pair is stabilized or not. This point 
has casted a serious doubt on the scenario, but in the 
course of the investigation on UPta , it has been pointed 
out that an odd-parity triplet pair is realized only when 
the inversion center exists external to the /-shell ions 
[l(ij |. This has been further examined by the detailed 
group-theoretical analysis 111 and the estimation of pair- 
ing potential with the use of band-calculation results [3| . 
However, the discussions have been in a phenomenolog- 
ical level and it is highly required that the triplet pair 
induced by the Hund's rule coupling should be investi- 
gated from the microscopic viewpoint. 

In this Letter, we attempt to gain an insight into triplet 
pairing induced by Hund's rule coupling. First we re- 
consider the symmetry argument on triplet pair in the 
weak-coupling limit. It is again fomid that the odd-parity 
triplet pair induced by the Hund's rule interaction occurs 
only for the non Bravais lattice. In order to confirm this 
result from the microscopic viewpoint, we carefully ana- 
lyze an orbital degenerate model on the two-dimensional 
(2D) square and honeycomb lattices. It is emphasized 
that the pair state is determined unambiguously by di- 
agonalizing the pair correlation function based on the 
concept of off-diagonal long-range order Then, we 

can visualize the pair wavefunction in an unbiased man- 
ner, clearly indicating that the odd-parity triplet pair in 
the FM phase occurs as an out-of-phase combination of 
local triplets induced by Hund's rule coupling in the 2D 
honeycomb lattice with two sites in the unit cell. 

Let us start our discussion on the Cooper-pair ampli- 
tude in the lattice with inversion symmetry. The pair 
operator composed of a couple of electrons is defined in 
the second quantized form as 



(1) 



where 



is an annihilation operator of an electron 



with spin a in the orbital a on the site /i included in the 
unit cell i. We consider the Fourier transform of Ci^aa as 



n,k 



(2) 




FIG. 1: (a) Two-dimensional square and (b) honeycomb lat- 
tices. Thin lines in (a) denotes 8-site cluster, while in (b), 
they indicate 8-site cluster composed of 4 unit cells including 
two sites, 1 and 2. Thick arrows define vectors ai and a2. 

where akna- is an annihilation operator for an electron 
with momentum k and spin a in the n-th band. Note that 
A^J^ is the matrix element of the unitary transformation 
between orbitals and bands. 

Now let us consider the Cooper pair of which the total 
momentum is equal to zero in the weak-coupling limit. 
Then, the Cooper-pair amplitude is given by 

{Kaa,i.f3r)^Yl e^"' """^ ^ A^^'^" (afc„.a_fe„.) . (3) 

Note that the Cooper pair is formed only by a couple 
of electrons on the same electronic band, since the pair 
composed of electrons in different bands is not stable in 
general in the weak-coupling limit. 

Here we introduce an inversion operator V, which acts 
generally on the electron operator as 

where Pa denotes the parity for the inversion of atomic 
orbital, depending on the angular momentum. In this 
paper, we assume the pairing in the orbitals with the 
same angular momentum, indicating Pa=Pp. 

For the Bravais lattice in which one site is included 
in the unit cell, typically the square lattice shown in 
Fig. 1(a), the inversion center is located on the site. 
Thus, we can simply suppress the indices fi and i> in 
this case. Combing Eqs. Q and Q, we obtain 

A^" = PaA-'^", (5) 

for the Bravais lattice. On the other hand, for the non 
Bravais lattice, /i is not equal to z/ in Eq. As a 

typical example, here we consider the honeycomb lattice 
composed of unit cells including two sites, as shown in 
Fig. 1(b). In this case, the inversion center is located at 
the center of two sites and we obtain 

At = Po.A-J;\ (6) 

First we consider the Cooper pair in the Bravais lattice. 
Note again that the index fi is suppressed in this case. 
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FIG. 2: Schematic views for wavefunctions of (a) odd-parity 
and (b) even-parity triplet pairs in the honeycomb lattice. 

After algebraic calculations, spin singlet (S) and triplet 
(T) pairs are, respectively, given by 

(•^ia, = E ^0^1'^ ■ - J«j)]A-'^"Af 5fc„, (7) 

and 

(•^i'a J/3) = E ' • - i2j)]A-'="A^"r,^„, (8) 

n,k 

where Skn = ((afenta-fcru) - {aknia-kn]))/V^, T^^ = 

»/V2, 

and Tj^^ — (aknia—kni) ■ As expected, spin singlet and 
triplet pairs have even- and odd-parity, respectively. 

Note here that the magnitude of on-site pair-amplitude 
vanishes for the odd-parity triplet pair, while for the sin- 
glet case, it depends on the symmetry of Cooper pair. As 
is well known, the Hund's rule coupling stabilizes the lo- 
cal triplet pair composed of electrons in different orbitals, 
inevitably leading to the orbital antisymmetric pair with 
finite on-site amplitude [l^ . However, such an even- 
parity triplet pair does not appear in the weak-coupling 
limit. In order to make this point clear, it is instruc- 
tive to consider the on-site pair-amplitude in the two- 
band model on the Bravais lattice. After simple algebraic 
calculations, we obtain (ciaaCifja) ^J2ki'^kiaa-k2a) for 
a^(3. Namely, the triplet pair induced by the Hund's 
rule coupling should be composed of electrons on dif- 
ferent bands, which is not stable in the weak-coupling 
limit We conclude that for the Bravais lattice, odd- 
parity triplet pair cannot be induced by the Hund's rule 
coupling. 

Next we turn our attention to the pair formation in the 
non Bravais lattice, typically the honeycomb lattice. For 
simplicity, we consider only the on-site pair amplitude. 
Using Eqs. Q and Q, for even-parity singlet and odd- 
parity triplet pairs, we obtain 

(<AL.ii/3) = (l/2)E(Ai2Ar;" + A-^"A^^)5;,„, (9) 
and 

(4a. n(3) = (1/2) E(Ai:A-/" - A-^"A^^)r,^„, (10) 

n,fc 

respectively. Note that in general, {(f>Jia ii/3) ^^^^ 
vanish. In order to understand this point, it is useful to 
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remark the relation 



^ilQ.il/g) 



[2/3/ 



(11) 



The pair satisfying this relation is composed of the local 
triplet pairs in the unit cell in an out-of-phase manner, as 
schematically shown in Fig. 2(a). Then, the odd-parity 
triplet pair induced by the Hund's rule coupling can ap- 
pear in the non Bravais lattice. 

We believe that the above symmetry argument is use- 
ful, but it mentions nothing about the stability of the 
odd-parity pair. For instance, the in-phase combination 
of local triplet pairs with even-parity [see Fig. 2(b)] seems 
to dominate the odd-parity pair at a first glance, but the 
above discussion does not conclude which pair is stabi- 
lized. In order to complete the discussion, we cannot 
avoid to carry out some explicit calculations in an ap- 
propriate model. However, if we apply mean-field ap- 
proximations on the model in which attractive interac- 
tions are introduced just by hand, we may lose essential 
points. Thus, in this paper, we adopt the orbital degen- 
erate Hubbard model, which is widely believed to be a 
standard model for strongly correlated electron systems. 
Furthermore, in order to obtain unbiased results, we re- 
sort to the numerical method such as exact diagonaliza- 
tion. For the purpose to analyze the paring symmetry 
in the numerical calculation, we apply the method of the 
optimization of pair correlation function 14] . 

Here we take the eg-orbital Hubbard model as 

H = - ^"/3(cLaaCj^/3<T + h.c.) 

(7,(7' .,a.^ (3 

+ •^'Zm^LaT'^LaiCiM/3iCiM/3T' (12) 

where denotes a pair of nearest-neighbor sites, 

V indicates a vector connecting the nearest-neighbor 
sites, a (/3) denotes x^—y'^ (3z^— r^) orbital, Pif^aa^ 

'^luaa^'^p-aa , and Piy,a^Y.a P'^t^aa ■ In thc first term, t^^ 
is a hopping amplitude depending on orbitals, hopping 
directions, and lattice type. Hereafter along the y 
direction (see Fig. 1) is taken as the energy unit. In 
the Coulomb interaction terms, C/, U' , J, and J' de- 
note intra-orbital, inter-orbital, Hund's rule, and pair- 
hopping interactions, respectively. Note the relations of 
U^U'+J+J' and J^J'. 

First let us briefly discuss the ground-state phase dia- 
grams for two types of 2D lattices with 8 sites (see Fig. 1). 
In Fig. 3(a), the phase diagram for the 8-site honeycomb 
lattice is shown for n=0.75. When J is increased, the 
FM phase appears, since it can gain the kinetic energy 
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FIG. 3: (a) Ground-state phase diagram for 8-site honeycomb 
lattice with n—0.75. Note that the region for J>U' is unphys- 
ical. (b) Maximum and second maximum eigen-values of P 
as a function of U' — J in the FM phase. CoefHcients of the 
odd- and even-parity triplet pair correlation functions for (c) 
a=(0,0), (d) a=ai, (e) a=a2, and (f) a=ai+a2, for U' — J—l. 



for large J. This is quite general in the orbital degenerate 
model. For n—1 (not shown here), the FM phase is also 
found, but the region becomes narrow: It does not touch 
the line of J=U'. For the 2D square lattice with 8 sites, 
the phase diagram for n—1.5 is quite similar to Fig. 3(a) 
(seeRef. [l^). For rt=l, the phase diagram is also similar 
to that of the honeycomb lattice with n=l. For n—0.75, 
the FM phase disappears in the square lattice. 

Now we focus on the pair in thc FM phase to clar- 
ify the symmetry of triplet pair induced by thc Hund's 
rule coupling. For the purpose, we measure the triplet 
pair correlation function, defined as P=($$t), where $ 
is a triplet pair operator given by a linear combination 
of (fifla ji//3- "1"^^ coefficient of each component includes 
information about pairing symmetry. The pair wavefunc- 
tion should be determined as the eigen function with the 
maximum eigenvalue of P |l4j . The pairing state thus 
determined, called "optimized pair" , is the most proba- 
ble candidate for the Cooper pair in actual systems. 

Before proceeding to the exhibition of numerical re- 
sults, let us briefly discuss the meaning of the optimized 
pair from the conceptual viewpoint of off-diagonal long- 
range order |l3l ]. In principle, the occurrence of super- 
conductivity is detected when the maximum eigenvalue 
Pmax of P becomes the order of N, where N is the num- 
ber of electrons. In a small cluster calculation, the possi- 
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ble superconducting pair state should be defined by the 
eigenstate with pmax, but in order to prove the existence 
of ofF-diagonal long-range order, it is necessary to show 
that Pmax actually becomes the order of N with increas- 
ing the system size. However, at some cluster size, a 
possible pairing state is determined without ambiguity 
by the eigenstate with pmax- 

In Fig. 3(b), we show the eigenvalues of P in the FM 
phase of the honeycomb lattice for n—Q.75. For small 
U' , the eigenstate of P corresponding to /Omax exhibits 
odd-parity, while the even-parity state appears for large 
value of U' 0. As discussed below, the symmetry of 
those pairing states are consistent with the schematic 
views shown in Figs. 2(a) and (b). In the fully spin- 
polarized phase, the relevant interaction is U'—J. Then, 
small U'—J corresponds to the "weak-coupling" region 
in the FM phase. In that sense, it seems natural that 
the odd-parity state appears for small U'—J, while the 
even- parity occurs for large U'—J. 

In order to visualize the paring symmetry, we explicitly 
express the pair operators as 

* = y^('<gla.lgCilaT'^i+al/3T + ^ la.2gCilaf Ci+a2/jT 
i,a 

+ 9'2a4/3Ci2QT'^i+al/3T + V2a,2f3'^i2a'\ Ci+a2P1 

+ 9'la,2aCilaTCi+a2QT + K/3,2/3Cil/3TCi+a2/3T): (13) 

where a is a vector connecting two unit cells. See 
Fig. 1(b) for the definitions of a in the honeycomb lattice. 

In Figs. 3(c)-(f), we depict the coefficients V'pa z^b 
each a and U'—J=l The symmetry argument in the 
weak-coupling limit has suggested that the odd-parity 
triplet pair should sat isfy the relations, (p^^ ip=~V2a 2f3 
and ^ia.2fj—~^2a.ii3 uM- Fig- 3(c), these relations 
are actually satisfied in the pair correlation inside the 
unit cell for a=(0, 0), indicating unbiased evidence for the 
odd-parity triplet pair. As shown in Figs. 3(d)-(f), the 
above relations for a=(0,0) are also observed for a^(0,0) 
and the pair amplitudes for a.y^(Ofl) are in the same or- 
der as those for a=(0,0). These results indicate that the 
triplet pair spatially extends to gain the kinetic energy 
by keeping the odd-parity symmetry. Here note an in- 
phase relation in the coefficients for a2 direction, while 
an out-of-phase relation for ai direction, suggesting a 
node structure between the unit cells along ai direction. 
Namely, the p-wave-like pair is appearing in this case. 
Note also another relation ffa,2a ~ ~^ii3 2/3 o*^*^" 
parity pair, indicating the node of the pair wavefunction 
in the orbital space. 

Finally, we briefly comment on actual compounds from 
the present result. First, the odd-parity triplet pair in 
the FM phase induced by the Hund's rule coupling is 
expected to occur in some /-electron compounds such as 
UGe2 and URhGe, although further quantitative calcula- 
tions are needed based on more realistic models. Second, 



the triplet pair in ruthenate may not be induced by the 
local Hund's rule interaction, since the ruthenate is de- 
scribed by the t2g-orbital degenerate Hubbard model on 
the 2D square lattice. Rather we may consider a scenario 
based on the pairing induced by the effective interaction 
among quasi-particles in the Fermi-liquid theory [20j . 

In summary, we have discussed microscopic aspects of 
odd-parity triplet pair. The weak-coupling analysis has 
suggested that the odd-parity triplet pair induced by the 
Hund's rule coupling does not appear in the Bravais lat- 
tice. We have carefully analyzed the triplet pair wave- 
function in the FM phase on the 2D square and honey- 
comb lattices. It has been clarified that the odd-parity 
triplet pair induced by the local Hund's coupling occurs 
only for the non Bravais lattice. 
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